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CHAPTER TWO 



ONE DIMENSIONAL WAVES IN SIMPLE SYSTEMS 
1 General solution to wave equation 



It is easy to verify by direct substitution that the most general solution of the one 
dimensional wave equation: 

can be solved by 



where /(^) and g{^) are arbitrary functions of ^. In the x, t (space,time) plane f{x — ct) 
is constant along the straight line x — ct — constant. Thus to the observer (x, t) who 
moves at the steady speed c along the positivwe x-axis., the function / is stationary. 
Thus to an observer moving from left to right at the speed c, the signal described 
initially by f{x) at t = remains unchanged in form as t increases, i.e., / propagates 
to the right at the speed c. Similarly g propagates to the left at the speed c. The lines 
X — ct ^constant and x -\- ct — contant are called the characteristic curves (lines) along 
which signals propagate. Note that another way of writing (1.2) is 



(t>{.x,t) = f{x 



ct) + g{x + ct) 



(1.2) 



(p{x, t) = Fit - x/c) + Git + x/c) 



(1.3) 



Let us illustrates an application. 



2 Branching of arteries 



References: Y C Fung : Biomechanics, Circulation. Springerl997 
M.J. Lighthill : Waves in Fluids, Cambridge 1978. 
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Recall the governing equations for pressure and velocity 



The two are related by the momentum equation 



= c^T^ (2.2) 



The general solutions are : 



Since 



and 



dp 



p = pj^{x — ct) + p-{x + ct) (2.4) 
u — Uj^{x — ct) + U-{x + ct) (2.5) 



dp / , / 



du- , , 

p-g^ = -pcu^ + pcu_ 



where primes indicated ordinary differentiation with repect to the argument. Equation 
(2.3) can be satisfied if 

p+ = pcu+, p_ = —pcu-. (2.6) 
Denote the discharge by Q = uA then 

ZQ± = u±A = ±p± (2.7) 

where 

Z^.g^^ (.8) 

is the ratio of pressure to flux rate and is call the impedance. It is the property of the 
tube. 

Now we examine the effects of branching; Referening to figure x, the parent tubes 
branches into two charaterized by wave speeds Ci and C2 and impedaces Zi and Z2. 
An incident wave approaching the junction will cause refiection in the same tube 

p = pi{t - xjc) + prit + xjc) (2.9) 
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and transmitted waves in the branches are pi{t — x / ci) and P2{t — x/c2). At the junction 
X = we expect the continuity of pressure and fluxes, hence 

Pi(t)+Pr(t)^Pl(t)^P2(t) (2.10) 

Define the rcfiection coefficient R to be the amphtude ratio of reflected wave to incident 
wave, then 

^ ^ ^ l~ + a) ,2 121 

Similarly the tranmission coefficients are 

Note that both coefficients are constants depending only on the impedances. Hence the 
transmitted waves are similar in form to the incident waves except smaller by the factor 
T. The total wave on the incidence side is however very different. 

3 Waves in an infinite domain due to initial distur- 
bances 

Recall the governing equation for one-dimensional waves in a taut string 

— -c— = 0, -oc<x<oo. (3.1) 
Let the initial transverse displacement and velocity be given along the entire string 

^x,0) = /(x) (3.2) 

^(x,t) = ^(x), (3.3) 

where ]{x) and g(x) are non-zero only in the finite domain of x. At infinities x ±oo, u 
and du/dt are zero for any finite t. These conditions are best displayed in the space-time 
diagram as shown in Figure 1. 
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ft AX 



7 7 

M=/fx) U=g(x) 



Figure 1: Summary of the initial-boundary- value problem 

In (3.1.1) the highest time derivative is of the second order and initial data are 
prescribed for u and du/dt. Initial conditions that specify all derivatives of all orders 
less than the highest in the differential equation are called the Cauchy initial conditions. 

Recall that the general solution is 

u = 0(0 + i^ir]) = (t>{x + ct) + i^ix - ct), (3.4) 

where 4> and ip are so far arbitrary functions of the characteristic variables ^ — x + ct 
and T] — X — ct respectively. 

Prom the initial conditions we get 

u{x,0) — 4>{x) + ip{x) — f{x) 
^{x,0) = c<f>'{x)-c,l^'{x)^g{x). (3.5) 

The last equation may be integrated with respect to x 



- f g(x')dx' + K, (3.6) 

C An 



where K is an arbitrary constant. Now (j) and tjj can be solved from (3.1.6) and (3.1.7) 
as functions of x, 

<t>{x)^\[f{x) + K] + ^J\{x')dx' 



i^{x)^\[f{x)-K]-^J\{x')dx\ 
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(x,t) 
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x-ct 



.4 



/ 

Domain of 
dependence 



c \ 



Figure 2: Domain of dependence and range of influence 



where K and xq are some constants. Replacing the arguments of by x + and of 
by X — and substituting the results in li, we get 



which is d'Alembert's solution to the homogeneous wave equation subject to general 
Cauchy initial conditions. 

To see the physical meaning, let us draw in the space-time diagram a triangle formed 
by two characteristic lines passing through the observer at x,t, as shown in Figure 2. 
The base of the triangle along the initial axis t — Q begins x — ct and ends at x-\-ct. 
The solution (3.1.9) depends on the initial displacement at just the two corners x — ct 
and x + ct, and on the initial velocity only along the segment from x — ct io x + ct. 
Nothing outside the triangle matters. Therefore, to the observer at x, t, the domain 
of dependence is the base of the characteristic triangle formed by two characteristics 
passing through x, t. On the other hand, the data at any point x on the initial line t — 
must influence all observers in the wedge formed by two characteristics drawn from x, 
into the region oi t > 0; this characteristic wedge is called the range of influence. 

Let us illustrate the physical effects of initial displacement and velocity separately. 





x+ct 



(3.7) 
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Figure 3: Waves due to initial displacement 

Case (i): Initial displacement only: f{x) ^ and g{x) = 0. The solution is 

u{x,t) = ^f{x - ct) + ^f{x + ct) 

and is shown for a simple f{x) in Figure 3 at successive time steps. Clearly, the initial 
disturbance is split into two equal waves propagating in opposite directions at the speed 
c. The outgoing waves preserve the initial profile, although their amplitudes are reduced 
by half. 

Case (ii): Initial velocity only: f{x) = 0, and g{x) ^ 0. Consider the simple example 
where 

g{x) = go when < b, and 
= when > 0. 

Referring to Figure 4, we divide the x ^ t diagram into six regions by the characteristics 
with B and C lying on the x axis at x — —b and +b, respectively. The solution in various 
regions is: 



in the wedge ABE; 
in the strip EBIF; 



1 go 
— J godx' ^ ^{x + ct + b) 



j-X+Ct 

u = — godx' = got 

Jx-Ct 
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in the triangle BCI; 



in the wedge FIG; 



Figure 4: Waves due to initial velocity 



2c 7_b c 



u 



— f go dx' = —{b — X + ct) 
2c Jx-ct 2c 



in the strip GICH; and 



in the wedge HCD. The spatial variation of -u is plotted for several instants in Figure 
4. Note that the wave fronts in both directions advance at the speed c. In contrast to 
Case (i), disturbance persists for all time in the region between the two fronts. 

4 Reflection from the flxed end of a string 

Let us use the d'Alembert solution to a problem in a half infinite domain a; > 0. Consider 

a long and taut string stretched from a; = to infinity. How do disturbances generated 
near the left end propagate as the result of initial displacement and velocity? 
At the left boundary x — Q must now add the condition 



u = Q, a; = 0, t > 0. 



(4.1) 



In the space-time diagram let us draw two characteristics passing through t. For 
an observer in the region x > ct, the characteristic triangle does not intersect the time 
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axis because t is still too small. The observer does not feel the presence of the fixed end 
at X = 0, hence the solution (3.1.9) for an infinitely long string applies, 

"X + Ct 



1 1 f^+'^^ 

^ = o [-^(^ + '^^) + -^(^ ~ '^^)] + ^ / 9{r)dr, X > ct. 



(4.2) 



But for X < ct, this result is no longer valid. To ensure that the boundary condition 
is satisfied we employ the idea of mirror reflection. Consider a fictitious extension of 
the string to — oo < x < 0. If on the side x < the initial data are imposed such that 
f{x) = — /(— x), g{x) = —g{—x), then u{0,t) = is assured by symmetry. We now 
have initial conditions stated over the entire x axis 

u{x,0) = F{x) and Ut{x,0) = G{x) — oo < x < oo, 

where 



F{x) 



f{x) if a; > 
-f{-x) ifa;<0 



g{x) if X > 
—g{—x) if a; < 0. 



G{x) = 

These conditions are summarized in Figure 5. Hence the solution ior < x < ct is 



u ^ I [F{x + ct) + F{x - ct)] + ^( r + [ ^ ] G{x')dx' 

^ \Jx-ct Jo / 

= I [/(^ + Ct) - f{ct - x)] + ^ (£ ^ + 9{x')dx' 

1 1 /"Ci+a; 

= ^ [/(^ + Ct) - f{ct -x)] + - g{x')dx'. (4.3) 

The domain of dependence is shown by the hatched segment on the x axis in Figure 6. 



5 Forced waves in an infinite domain 



Consider the inhomogeneous wave equation 
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M=0 



u =c u 

tt XX 



u^-f(-x) 
u ^-g(-x) 



u^fix) 
u=g(x) 



Figure 5: Initial-boundary- value problem and the mirror reflection 



X+Ct=XA-Ct„ . 

^ °^ t 



X-Ct=X^-Ctg 



Figure 6: Reflection from a fixed end 
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where h{x, t) represents forcing. Because of linearity, we can treat the effects of initial 
data separately. Let us therefore focus attention only to the effects of persistent forcing 
and let the initial data be zero, 

du 



u{x,0) = 0, 
The boundary conditions are 



dt 



= 0, (5.1) 

t=o 



— > 0, \x\ — > oo. (5.2) 
Let the Fourier transform of any function f{x) and its inverse f{a) be defined by 

/oo 
f{x)e-"^^dx (5.3) 
■oo 



1 



oo 



/(a:) = — y f{a)e'''^da (5.4) 

The transformed wave equation is now an ordinary differential equation for the transform 
of u{x, t), i.e., u{a, t), 

d Uj ey ey — 

— - + c^'a^'u ^ h t>0 

at'' 

where h{a, t) denotes the transform of the forcing function. The initial conditions for u 
are: 

tZ(a, 0) = /» = 0, ^^^^-g{a)^0 

Let us hide the parametric dependence on a for the time being. The general solution 
to the the inhomogeneous second-order ordinary differential equation is 

U = CiUiit) + C2U2{t) + ^ ^ [ui{T)u2{t) - U2{T)u2{t)\ dr, (5.5) 

where Ui and U2 are the homogeneous solutions 

= e-^"^* U2 = e''^"' 

and W is the Wronskian 

W — U1U2 — U2u'i — 2iac — constant. 
The two initial conditions require that Ci = C2 = 0, hence the Fourier transform is 



u 
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The inverse transform is 

"t 1 POO 



U{x,t) = — [ dr- r ^'^^ [e^-Mt-r)) _ ^ia(x-c(t-r))^ (5 7) 

271 Jq 2c J_oo ioi 



The integrand can be written as an integral: 

1 1 rx+c(t—T) 



u = — 
27r 



/"t -j^ /"CXD PX+C^t — T) 

/ dr— / dah{a,T) / d^e*"^ 

Jo 2c Jx-c{t-T) 



By interchanging the order of integration 

■j^ pt r-x+c{t-T) 
' X — c{t—T) 

and using the definition of Fourier inversion, we get 



2c 



rt /•x+c{t—T) -I /"OO 

dr di— / dQ;e^"«/i(a,r) 

Jo Jx-c(t-T) 27r 



ti(x,i) = - / dr / deMe,T) (5.8) 

Jq Jx-c{t-T) 

The right-hand side is the integration of h over the characteristic triangle defined by 
the two characteristics passing through {x,t). in the x — t plane. Thus the observer is 
affected only by the forcing inside the characteristic triangle. 

For non-zero initial data u{x, 0) = f{x) and Ut{x, 0) = g{x), we get by linear super- 
position the full solution of D'Alambert 

<x,t) = ^[f(x + ct) + f(x-ct)] + -J ^d^giO 

1 pt l'X+c{t-T) 

+ dr dCHCr), (5.9) 

Jo Jx-C{t-T) 

The domain of dependence is entirely within the characteristic triangle. 



6 String embeded in an elastic surounding 

Reference : Graff : Waves in Elastic Solids 

If the lateral motion of the string is restrained by elastic springs along the entire 
length, the governing equation can be found from §1.1 by replacing the external pressure 
by the elastic restoring force —KV per unit length. 
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which can be written as 

1 d^V 

where 

Co 

6.1 Monochromatic waves 

For any hnearized wave problem, if the range of the spatial corrdinate x is (—00,00), 
and all coefficients are independent of x, then the first task is to examine the physics 
of sinusoidal wave train of the form: 

where A — \A\e'^'^^ is a complex number with magnitude \A\ and phase angle 0^- After 
examining the physical meaning of this special type of waves, it is possible to use the 
principle of superposition to construct more general solutions. It is customary to omit 
the symbol IR.— "the real part of for the sake of brevity, i.e., 

V{x, t) = ^6^'=^-^'"* (6.5) 

First of all a few definitions about sinusoidal waves in general. We shall call 

e{x,t) = kx - uut (6.6) 

the wave phase. Clearly the trigonometric function is periodic in phase with the period 
271. In the X, t plane, V has a constant value along a line of contant phase. In particular, 
9 = 2n7i, (n = 0, 1, 2, ■ ■ ■ ) correspond to the wave crests where V = \A\ is the greatest. 
On the other hand, 9 — {2n + l)n, (n = 0, 1, 2, • • • ) correspond to the wave troughs 
where V — —\A\ is the smallest. \A\ is half of the separation between adjacent crests 
and troughs and is called the wave amplitude; we also call A the complex amplitude. 
Clearly || represents the number of phase lines per unit distance, i.e., the density of 
phase lines, at a given instant; it is called the wavenumber, 

wavenumber — k — — 

ox 



dx^ 



-V 



(6.2) 



(6.3) 



(6.7) 
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On the other hand — 1| represensts the number of phase hnes passing across a fixed x 
per unit time; it is called the wave frequency. 



wave frequency —co- 



de 
di 



(6i 



To stay with a particular line of contant phase, say a crest, one must have 

d9 — kdx — ujdt — 
namely one must move at the phase velocity, 



dx 


CO 


c—-r 
dt 


6i=constant ^ 



(6.9) 



Now back to the string. Substituting (6.10) in (6.1) we get 

,2~ 



,2 ^ 

1 ci 



V = Q 



or 



uj = Co\ A;2 + 



K 
T 



The phase speed is 



W\ 



K 



1/2 



> Cn 



(6.10) 



(6.11) 



(6.12) 



See figure 6.1. Note that, due the stiffening by the lateral support, the phase speed 
is always greater than Co, and decreases monotonically with the wave number. Longer 
waves (small \k\) are faster, while shorter waves (larger |A;|) are slower. As increases, 
c approaches the finite limit Cq for the shortest waves. In general a sinusoidal wave whose 
phase velocity depends on the wavelength, i.e., a; is a nonlinear function of A;, is called 
a dispersive wave, and (6.11) or its equivalent (6.20) is called the dispersion relation. 

An interesting physical feature for dispersive waves in general can be found by su- 
perposing two trains of sighltly different frequencies and wave numbers: 



where 



V = A ^g*('="'"^-'^"''*) _|_ /le**-*^ 



k^^k±k', u^^uj{k^), k' <^k. 



(6.13) 



(6.14) 
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Figure 7: Phase and group velocities of a string in elastic surrounding. Point of Sta- 
tionary phase ko. 



Let us approximate o;^ by 



= u(k) ± k''^ + 0(k'^) 
dk 



then 

Y — j^^ikx-iujt ^-ik'{x-{dw/dk)t) _|_ ^-ik'{x-{dw/dk)t)^ _ j^^^ ^^^ikx-iuit 

where 

2cosk' {x - {du;/dk)t) (6.16) 

The factor exp{ikx — iuit) is called the carrier wave and A{x,t) the envelope. Thus the 
result is a sinusoidal wave train with a slowly varying envelope which has a very long 
wavelength 27r/A;' ^ 27r/A; and moves at the group velocity 



did 

group velocity = c„ = — (6.17) 

dk 



which is in general different from the phase velocity for dispersive waves. 

In our string problem, the group velocity is easily found from the dispersion relation 
(6.11) 



,2 



c, = —J^^^= = ^^—, hence ^ = %<1 (6.18) 
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Thus the group velocity is always smaller than the phase velocity, and increases with 
the wavenumber from for the longest wave to the finite limit Cq (equal to the phase 
velocity) for the shortest waves. 



When u) > u)c where 

Ik 

T' 

is called the cutoff frequency ujc, k is real 



(6.19) 



k^±\k\, with l^l=(^^-^J (6-20) 

being the real and positive square root. V{x., t) is a propagating wave, with the plus 
(minus) sign corresponding to right- (left-) going wave. If a; < cUc, k — ±iK is imaginary, 
with K being the postive root: 

,2\ 1/2 



Then e^'^^ — e^'^^. For boundedness one must choose the minus (plus) sign for x > 
{x < 0). Oscillations are localized or evanecent; there is no wave radiation. 

As a simple application, (6.4) is the response in a semi-infinite string forced to 
oscillate at the left end x = 0, 

V{0,t) = ^{Ae-''^'}. (6.22) 

li cu > CUc, then 

V{x, i) = 3? [Ae'^k\x-i^t^ (g_23) 

where \k\ is defined by (6.20). The requirement that waves due to a local disturbance 
can only radiate outwards is called the radiation condition. More will be said on it in 
Appendix B. 

If a; < a;c, we must require boundedness at infinity so that 

V{x,t) = ^Ae-''''-''^\ x>0 (6.24) 

At the cutoff frequency. A; = 0, is constant in x; the infinitely long string would 
oscillate in unison. This absurd result signfies the breakdown of the the linearized 
theory. 



2.6. STRING EMBEDDED IN AN ELASTIC SUROUNDING 



17 



6.2 Energy transport 

Along a unit length of the string the kinetic energy density is, 

The potential energy in any segment dx of the string is the work needed to deform it 
from static equilibrium. The part due to lengthening of the string against tension is 



T{ds - dx) = T\ll+ (^— j dx -Tdx^- (^— j dx 



Adding the part against the springs, the total potential energy per unit length is 

We now calculate their time averages. If two time-harmonic functions are written in the 
complex form: 

a = 3? (Ae-^'^*) , 6 = 3? {Be''''*) (6.27) 

the time average of their product is given by 

ft+— 

^ = — / ahdt= -lk{AB*) = -lk{A*B) (6.28) 
27r 2 2 

Using this recipe, the period-averaged energy densities are. 



PE^^^ {ikAe^f + ^{^Ae^f = {^e + 1^1 

where 9 = \k\x — cut is the wave phase. Since 



2 T f,^ K 



\rho V T 
we have 

P 



PE=^uj^\A\^ 
4 

Thus the total energy is equally partitioned between kinetic an potential energies, and 

P. .21 Al2 



E = KE + PE=-uj^\AY (6.29) 



2.7. DISPERSION FROM A LOCALIZED DISTURBANCE 



18 



Now the averaged rate of energy influx across any station x can be calculated from 
the product of transverse component of the tension and the transverse velocity, 

This power input is transported from left to right by the wave. Hence the speed of 
energy transport is the group velocty. This result is quite general for many physical 
problems and is not limited to the springs. 



7 Dispersion from a localized initial disturbance 

The solution for monochromatic waves already shows that waves of different wavelengths 
move at different velocities. What then is the consequence of an initial disturbance? 
Since a general initial disturbance bounded in space can be represented by a Fourier 
integral which amounts to the sum of infinitely many sinusoids with a wide spectrum, 
we shall employ the tools of Fourier transform. 
In addition to the governing equation: 

1 d^V d^V KV 



c2 df^ dx^ T ' 
we add the initial (Cauchy) conditions 



oo<a;<oo, t > (7.1) 



V(.,0)=/W, (7.2) 

Let us define the Fourier transform and its inverse by 

/oo 1 poo 

e-''^fix)dx, fix) = — e^'^mdk (7.3) 



The transforms of (7.1) is 



c2 dt^ T ' 

and of the intitial conditions are. 



-k^V-—, t>0 (7.4) 



V{KQ) = m, ^^^ = (7.5) 
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The solution for the tranform is 

V{k,t) ^ f{k) cos ut (7.6) 

where 



The Fourier inversion is 



u; = Co^k^ + ^ (7.7) 
1 f°° 

V{x,t) = — dkf{k)e'^'' cos ujt (7.8) 

Any real function f{x) can be expressed as the sum of an even and an odd function 
of X. For simplicity let us assume that f{x) is even in x so that f{k) is real and even in 
k, then 

1 /"^ 

V{x.,t) = — I dkf {k) cos kx cos ujt 
which can be manipulated to 

V{x, t) = — 3? / dkf{k) [e'^^-i^' + 6*'=^+^^*) dk (7.9) 
27r Jo 

The first term in the integand represents the right-going wave while the second, left- 
going. Each part correponds to a superposition of sinusoidal wave trains over the entire 
range of wave numbers, within the small range {k, k + dk) the amplitude is f{k)/2. The 
function f{k)/2 is called the Fourier amplitude spectrum. In general explicit evalua- 
tion of the Fourier integrals is not feasible. We shall therefore only seek approximate 
information. The method of stationary phase is particularly useful here. It aims at the 
asymptotic approximation of the integral 

i{t) = I F{ky*^^^Uk (7.10) 

J a 

for large t. Let us first give a quick derivation of the mathematical result. Assume that 
F{k),(t){k) are ordinary functions of k. If t is large, then as k increases along the path 
of integration both the real and imaginary parts of the exponential function 

cos{t(f){k)) + i sm{t(f){k)) 

oscillates rapidly between -1 to +1, unless there is a point of stationary phase kg within 
(a,b) so that 

^^=0'(fc„) = O, a<ko<b. (7.11) 
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Figure 8: Neighborhood of Stationary Phase 

Then important contribution to the Fourier integral comes only from the neighborhood 
of ko- Near the point of stationary phase, we approximate the phase by 

and the integral by 

I{t) ^ F(fc,)e^*^(^°) ^'exp (^-^{k - kofcj)"{ko)^ dk 

With an error of 0{l/t), we also replace the limits of the last integral by ±oo; the 

justification is omitted here. Now it is known that 

oo r~ 

±^tk^'dk = J-e^'^" 



1/2 /^X 

+ ( - j , if /Co G (a, 6), 

(7.12) 

where the sign is + (or — ) if (j)"{ko) is positive (or negative). It can be shown that if 
there is no stationary point in the range(a,b), then the integral I{t) is small 

I{t) = o(^Y if ko^{a,b). (7.13) 



e 

V ^ 

It follows that 

I{t) = / F(A;)e**^('^) dk ^ F{ko)e'">'^^°^^'^'^ 
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Let us apply this result to the right-going wave 

1 

V+{x,t) = —^J dkf{ky^^'''''^-'''^dk (7.14) 



For a fixed x/t ^constant, we have 



K 

(f){k) = kx/t — u){k), with x/t = ^, u! = Cod k'^ + — (7-15) 

For an observer travelling at the speed x/t less than the maximum uj'{k), i.e., Cq, there 
is a stationary point ko at the root of 

x/t = J{k.) = ^ = (7.16) 

The root fco increases from zero for long waves to the maximum Cq for very short waves. 
If the observer is faster than Co, there is no stationary point. 
Since 

(t>"{ko) = -uj"{ko) < 0, (7.17) 

the final result is that 



The transform f{k) depends on the specific profile of the initial disturbance. For exam- 
ple, if 

which has an area S and characteristic width &, the Fourier transform is 

f{k) = Se-\''\^ (7.20) 

To a specific observer identified by the speed of travel x/t < c^, the approximate 
result can be viewed as a simple harmonic wave train with wave number koi frequency 
uj{ko) and amplitude 

A-^^^-^ I 2" (7 211 
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Figure 9: Snapshots of wave field on the side x > 0. 

Since uj'ih^j increases from with increasing h to the finite maximum T / pCo = Cq, an 
observer faster than Cq sees no waves. However any observer slower than Cq is accom- 
panied by a train of progressive sinusoidal waves. The local wavelength ko is such that 
the group velocty matches the observer's speed. The faster the observer, the shorter 
the waves. If a snapshot is taken, then the shortest wave, whose phase velocity is the 
lowest are seen at the front which moves as fast as the crests of the shortest waves. The 
crests of longer waves advance much faster than the local envelope. Because / (k) is the 
greatest at A; = 0, the envelope of the longest waves which stays near the source, is the 
biggest. The envelope of the shortest waves is lower in amplitude and spreads out with 
the wave front. The entire disturbance attenuates in time as See figure 7 for an 

overview. 

Note also that very near the wave front, Cg — > Cg', the second derivative (t)"{ko) — 
—uj"{ko) vanishes. Hence the asymptotic formula breaks down. A better approximation 
is needed, but is omitted here. (See C. C. Mei, 1989, Applied Dynamics of Ocean Surface 
Waves) . Finally we examine the propagation of wave energy in this transient problem. 
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Using (7.21) the local energy density is: 

At any given t, the waves between two observers moving at slightly different speeds, 
Cg{ki) and Cg{k2), i.e., between two points Xi/t = Cg{ki) and X2/t = Cg{k2) are essentially 
simple harmonic so that the total energy is 



/ dxE^ I 

J X1 J X 



Ant u"{ko) 
uj"{ko)dko 



' X\ J Xl 

Since x — uj'{ko)t for fixed we have 

dx 

T 

Now for X2 > xi, k2 > ki, it follows that 

/ dxE — / dko— — (/(/co))^ = constant 
Jxr Jki 47r 

Therefore the total energy between two observers moving at the local group velocity 
remains the same for all time. In other words, waves are transported by the local group 
velocity even in transient dispersion. 



8 Scattering of sinusoidal waves 

If along a long rod there are some localized inhomogeneities, an incoming train of sinu- 
soidal waves will be partly reflected and partly transmitted. The scattered signals tell 
us something about the scatterer. To determine the scattering properties for a known 
scatterer is called the scattering problem. To determine the scatterer from the scattering 
data is called the inverse scattering problem. We shall only study the former. 

Various mathematical techniques are needed for different cases: (i) Weak scatterers 
characterized by small amplitude relative to the wavelength, or slow variation within a 
wavelength, (ii) Strong scatterers if their dimensions are comparable to the wave length. 

8.1 Weak scattering 

Let the long rod have a slightly nonuniform cross section, 

S{x) = So{l + ea{x)), e < 1. (8.1) 
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where a{x) diminishes to zero at x ~ ±00. The wave equation reads 
For sinusoidal waves 

u{x,t)^U[U{x)e-'^^] 
the spatial factor is governed by 

^s'^)+pu;''SU^O (8.3) 



dx \ dx J 

Then the solution can be sought by the perturbation method 

U{x) = Uo{x) + eUi{x) + e^U2{x) + ■■■ (8.4) 

Substituting this into the governing equation 

ESo-^ |(1 + ea(x))^ [Uo{x) + eUi{x) + ^U2{x) + • • •] | 

+puj^So{l + ea{x)) [Uo{x) + eU^{x) + e^U2{x) + •••]= (8.5) 

and equating the coefficients of like powers of e to zero, we get from the terms of order 

0(eO): 

ESo'^ + puj''SoU^ = Q (8.6) 
The solution is simply the incident wave 



C/o = Ae''=^, with k^u^fp|E (8.7) 



At the order 0(e) we get 



~dx^ ^^^^ \ ~dx^ J dx~dx ^ 

or, 

^ + = (8.8) 

Thus the scattered wave is at most of the order e and is governed by an inhomogeneous 
equation. Let us define the fundamental solution (Green's function) by 

f^li^ + k^Gix, x') = Six - x') (8.9) 
dx^ 
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G ~ outgoing waves at ± oo 
It will be shown in Appendix A that 

Qik\x—x'\ 



G{x,x') 



2ik 



].10) 



Ul) 



By straightforward differentiation, it can be verified that 

da{x') dUo{x') 



Ui(x) = - 



ax' ax' 2ik ./_^ ax ax 



±l^ik\x-x'\^^> (8.12) 



Since Uo{x') — e**^^', we have 

A r da{x') 



2 ./_oo dx' 



A 



Akx 



da 
dx' 



dx' + e 



' — oo 

A ^ 



f 



°° dn 
ikx I ^^^ik^'dx' 

dx' 



(8.13) 



2 
A 



e*''^a(x) + e 



-ikx 



d^ ^2ikx' j^f 

dx' 



ikx I „—ikx 



a{x)e^'''' + e 



A 
dx 



e""""' dx 

- [aixy^^'^'^ dx' - a2ike''^''' 



dx' 



= ikAe-''''' / a{x')e^'''='' dx' 

J X 

Far to the right we get 



(8.14) 



t/i = 0, x^oo. (8.15) 

Thus on the transmission side modification of the incident waves can at most be of the 
order O(e^). Far to the left x ~ — oo, 



U^eUi^ Ae-'^'' (^iek J a{x')e'''''''' dx'^ , x 



-oo. 



Let us define the reflection coefficient by 

/oo 
a{x')e^'''''' dx' , 
■oo 

For instance if a is a hump of unit total dimensionless area and length b, 

aix] — 



7r(x2 + 62) 



then 



R — iekbe 



-2kb 



(8.16) 



(8.17) 



(8.18) 



(8.19) 
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Figure 10: Reflection coefficient given by ( 8.19). 



Hence reflection is small if kb is large, i.e., long and gentle obstacles are not effective 
in reflecting short waves . For very long waves or a short obstacle, reflection is also 
vanishingly small. The maximum reflection coefficient is R — ieao/2 when kb — 1/2. 
Note that far to the right, C/i(x) — > as x — > oo. The transmission vwave is 



(8.20) 



The tranmission coefficient is T = 1 + O(e^). This is consistent with the general law of 
energy conservation + |Tp = 1 to be proven later. 

Homework: Check if the governing equation at O(e^) for U2 is: 



+ puj'SqU2 = ESo 



da dUo da dUi 



dx dx dx dx 



.21) 



Carry out the solution to find the transmission coefficient T from the amplitude of 
Uq + eUi + at x ~ 00, and show that energy is conserved to the order O(e^). 
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8.2 Strong scattering of long water-waves by a shelf 

Consider an ocean bottom with a step- wise variation of depth. 

hi, X < —a; 

h = { h2, -a<x<a; (8.22) 
hs — hi, X > a 

How does the step interupt the propagation of an incident wave of unit amphtude 
arriving from x ~ — oo? 

In each zone of contant depth {i = 1, 2, 3), the shaUow water equations read: 

""•+9?? = (8.24) 



For sinusoidal waves 



we get 



dt dx 



Ci-Vie"^\ Ui^Uie-"'' (8.25) 



-iujr]i + h~^ = (8.26) 



dx 

-iuUi + g^^O (8.27) 



or 



72 

"^"^^ + k^rii = 0, where ki = ^= (8.28) 



dx^ ' " ' ' 
At a junction, the pressure and the flux must be equal, hence 

rji = 772, and hi^ = /i2-r^, x = -a; (8.29) 
ax ax 

?72 = ?73, and h2— = hi — , x ^ a. (8.30) 
The forms of the solutions in each zone of constant depth are: 

7]i = e^^'i(^+«) + i?e-^*^i(^+'^\ X < -a; (8.31) 
7^2 = Ae^*^2x ^ 5e-ife2x^ -a<x<a (8.32) 
= re**=i(^-«), X > a (8.33) 
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The reflection and transmission coefficients R and T as well as A and B are yet unknown. 
Applying the matching conditions at the left junction, we get two relations 

1 + R = Ae-'''^" + Be'^^" 

hhil -R) = k2h2{Ae-'^^'' - Be'^^"). 
Similarly the matching conditions at a; = a gives 

k2h2{Ae'^'"' - Be-'''^^) = hh^T. 
These four equations can be solved to give 

4s 



R 



_j_ g^2g— 2jfc2a ^^2g2i/c2a 

''^ _j_ g^2g— 2jfc2a ij^2g2ifc2a 

A = ^e-'^^" (1 + s) 

B = -e'^"^" (1 - s) 
2 ^ ' 

where 

kihi hi c\ 



^.2^2 V ^2 C2 

The energy associated with the tranmsitted and reflected waves are : 

T 



4s2 + (1 -s2)2sin2 2fc2a 
(1 - s^) sin^ 2k2a 



4s2 + (1 - s2)2 sin^ 2A;2a 
It is evident that \R\^ + \T\^ = 1. 

Over the shelf the free surface is given by 

2s [(1 + s)e^''2(^^-a) + (1 _ s)e-»fc2(^-a)j 

^2 = (1 + s)2e-ifc2a _ (1 _ s)2eifc2a (S'^^) 

Recalling the time factor e"*"^*, we see that the free surface over the shelf consists of two 
wave trains advancing in oppopsite directions. Therefore along the shelf the two waves 
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can interfere each other constructively, with the crests of one coinciding with the crests 
of the other at the same moment. At other places the interference is destructive, with 
the crests of one wave train coinciding with the troughs of the other. The envelope of 
energy on the shelf is given by 

2 [cos^ k2{x — a) + sin^ k2{x — a)] 

l'^' " 4s2 + (l-s)2sin2 2A;2a ^^"^^^ 

At the downwave edge of the shelf, x — a, the envelope is 

Note that the reflection and transmission coefficients are oscillatory in k2a. In par- 
ticulatr for 2A;2a = nn^n — 1,2,3..., that is 4a/A = n, \R\ — and \T\ — 1. The 
shelf is transparent to the incident waves, correspoinidng to the most constructive inter- 
ference and the strongest transmission Mininum transmision and maximum reflection 
occur when 2/c2a = (n — l/2)7r, or 4a/A = n — 1/2, when the interference is the most 
destructive. The corresponding transimssion and reflection coefficients are 



min|T|2 = 7^;—^' max|i?|^ = ^ (8.48) 



4g^ I _ (1 ~ 



See Figure 11 

The features of interference can be explained physically. The incident wave train 
consists of periodic crests and troughs. When one of the crests first strikes the left 
edge at X = a, part of the it is transmitted onto the shelf and part is reflected towards 
X ~ — oo. After reaching the right edge at a; = a, the tranmitted crest has a part 
reflected to the left and re-reflected by the edge x = —a to the right again. When the 
re-reflected crest arrives at the right edge the second time, its total travel distance is 
4a. If 4a is an integral multple of the wave length A2, the crest is in phase with all 
the other crests entering the shelf either before (the third, fourth, fifth, ... time) or 
after. Thus all the crests reinforce one another at the right edge. This is constructive 
interference, leading to the strongest tranmission to the right x ~ 00. On the other 
hand if 2A';2a = (n — l/2)7r or 4a/A = n — 1/2, some crests will be in opposite phase to 
some other crests, leading to the most destructive interference at the right edge, and 
smallest transmission. 
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Figure 11: Transmission and reflection coefficients for a rectangular shelf. Prom Mei, 
C.C., Applied Dynamics of Ocean Surface Waves. 

Prom the results in the following exercise, the reffected wave suffers a reversal of 
phase at the discontinuity (a crest turns into a trough and vice versa). Hence at the 
reflection edge a wave crest tranmitted onto the shelf returns from trough. 
It is then transmitted to the left of x — —a as trouigh. On the other hand, a crest that 
is reflected at x — —a appears as a crest. Therefore the interference is destructive. The 
resulting R, which is the sum of multiple reflections, is small. 

Homework: Scattering by a depth discontinuity: Consider the scattering of long 
water waves by a depth discontinuity at x — 0, i.e., h — hi,—oo < x < and h — 
h2,0 < X < oo. Show that the tranmssion and reflection coefficients are: 

i + (V^i)V2' ^ i + (V^i)V2 ^^-^^^ 

Thus the reflected wave suffers a reversal of phase when /i2/^i > 1, i-e, approaching the 
discontinuity from the shallow side. 
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9 General identities in scattering of monochromatic 
waves by arbitrary inhomogeneities 

Scattering due to inhomogeneities, caused by nonuniformities either in geometry or in 
material properties, requires the solution of ordinary differential equations with variable 
coefficients. In general one must resort to numerical means. For the sake of checking 
numerical accuracy and for gaining physical insight, identities which must be true are 
useful. They are often deduced by arguments typical in the derivation or the use of 
Green's theorem. 

We illustrate these identities by the example of an infinitely long rod with variable 
cross section S{x), 

dx {^^di) " ^~eW ^^-^^ 

E and p are taken to be constants. S{x) is nonuniform only in a finite neighborhood 
around a; = 0. Elsewhere S = So =constant. Consider monochromatic waves 

u(x,t) ^U(x)e-'''^ (9.2) 

so that U {x) is governed by the ordinary differential equation 

{SU')' + k^SU = 0, = ^ (9.3) 

E 

The boundary conditions depend on the source of the incident waves. 

Let us consider two problems: one with the incident wave from x ~ — oo; the other 
with the incident wave from x ~ oo. 

If the incident waves arrive from x ^ —oo, the asymptotic form of the solution must 

be 

t/i(x)-^_e''=^ + 5_e-''=^ = ^_(e^'=^ + i?ie-^'=^), kx ^ -oo (9.4) 

and 

Ui{x) ~ A+e''^^ = A.Tie''''' kx ^ oo (9.5) 

where Ti and Ri are the transmission and reflection coefficients which arc a part of the 
unknown solution of the left-incidence problem. Let us define the Jost function by 
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then 



/i(a;)~^e*'=^ + ^e-^'=^ fcx « -1 (9.7) 



and 

/i(x)~e^*^^ A;x>l (9.8) 

On the other hand, if the incident waves arrive from x ^ oo, the asymptotic form of 
the solution must be 

U2{x) ~ ft.e-^'^^ = ^)+T2e-^^^ kx ~ -oo (9.9) 

and 

U2{x) ~ a+e^^^ + fc+e-^'^^ = 6+ {R2e'^'' + e'^'^^) , fcx ~ oo (9.10) 

where T2 and i?2 are the unknown transmission and reflection coefficients of the right- 
incidence problem. Similarly we define another Jost function 

h^T-^ (9-11) 

so that it behaves as 

/2(x) - e-^*^^ kx^-00 (9.12) 

and 

f2ix)^^e-''^ + ^e^'^, kxr. oc (9.13) 

-1-2 J-2 

Since both /i and /2 satisfy (9.3), we have 

(5/0' + fc'^/i=0, and {Sf^' + k^Sh = Q, (9.14) 

Multiplying the first by /2 and the second by /i, then taking the difference, we get after 
partial integration, 

implying the following Wronskian identity 

W{h, h) = (/i/^ - /2/O = (9.15) 
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The constant C can be related to the asymptotic hmits at k\x\ » 1 where S ^ Sq. Far 
on the incident side, the Wronskian is 

(/1/2 ~ f2f[)x 00 

^ ikx I -ikx\ ( ■u\-ikx ^-ikx ( Akx '^^■^^ „-ikx 



zkx ^ ^ -zkx t^j^x -tkx _ -%kx %kx i 

Far on the transmission side, the Wronskian is 

(/1/2 ~ /2/l)a:~oo 

_ 2ik 

Since the two Wronskians are equal, we conclude that 

Ti = T2 (9.16) 

Thus, the complex transmission coefficients are equal for incident waves coming from 
either side, no matter how asymmetric the geometry may be in the near field! 

Next we consider the left-incidence problem and omitting the subscript " 1" . Multi- 
plying the first of (9.3) by the complex conjugate of U, i.e., U*, and taking the difference 
of the result with its complex conjugate, we get 

us{u*'y - u*{suy = ^ {s{uu*' - u*u')) = o (9.17) 

hence 

S{UU*' - U*U') = constant (9.18) 
By using the asymptotic expression (9.4) on the incidence side, we get 

[S{UU*' -SU*U')]^^_^ 

= So (A.e^'^^ + S.e-^'^^) {-ikA*_e-''''' + ikB*_^^^) 
-So {A*_e-'^'' + B*_(^^'') {ikA.e'^'' - ikB^e''^'') 
= -2ik\A_\^ + 2ik\B_? 
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Similarly by using the asymptotic expression (9.5) on the transmission side, we get 

[S{UU*' -U*U')]^^^ 

= So (A+e*'^^) {{-ik)A\e-''''') - So (A;e^*'=^) ((iA;)A+e^*^^) 
^ 2z A^iS'q I I 

Equating the two limits we conclude that 

= + (9.19) 

or 

+ |T^|2 ^ 1. (9.20) 

This is merely a statement of conservation flux : the flux rate of total scattered energy 
(reflected and transmitted) is equal to the flux rate of the incident wave energy; the 
speed of energy tranport being the same on both sides of the scatterer. Clearly the 
same energy conservation must hold for the right-incidence problem. 

|i?2|' + |T2|2 = 1 (9.21) 

If the scattering problem is solved numerically by, say, flnitc elements, it is necesssary 
that the computed scattering coefficients R and T satisfy the identities (9.16) (9.20) and 
(9.21). 

Homework: Extend the identities (9.16) and (9.20) when S{x) Si as kx ^ — oo 
and S{x) — > S2 as kx — > 00, where and S2 are two different constants. 

10 Refraction in a slowly varying medium 

Let us consider time-harmonic long waves in shallow water waves. The governing equa- 
tion is 

d / , dri\ cu"^ , , 

Consider a sea depth which varies slowly within a wavelength, i.e.. 
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Earlier analysis suggests that reflection is negligibly small. Thus the solution is expected 
to be a locally progressive wave with both the wavenunmber and amplitude varying much 
more slowly than the wave phase in x . Hence we try the solution 

7] = A(a;)e'^(^) (10.3) 

where 9{x) — cut is the phase function and 

JO 

k{x) = ^ (10.4) 
ax 

is the local wave number and A is the complex amplitude. Note the spatial rate of 
variation of the phase function, i.e., the wave number is in gneral not small, therefore 9 
itself is not a slowly varying function of x. Let us calculate the flrst derivative: 

Iaj<uU \ Hj<uU I 

and assume 

dA 

= 0{kL)-^ < 1 

In fact we shall assume each derivative of /i, A or A; is times smaller than /c/i, kA or 
k'^. Puthurmore, 



d ( drj\ ijj 



,2 



, , „ , , , , , dA\ d ( dA\ d(khA) u^A ,„ ^ 

, yh-T^]+—V= ik(ikh + h—)+ — ( h— ]+i \ V e'' ^ 

ax \ ax J g [ \ ax J ax \ ax J ax g 

Note that the complex amplitude A can be written as 

A = \A{x)\e'^^^''^ (10.5) 

The phase of A, 6a, is a slowly varying function of x. It can be considered as a part of 
the wave phase, although its spatial gradient is much smaller than k. 
Now let's expand 

A^ Ao + Ai + A2 + --- (10.6) 
with Ai/Aq — 0(/i), = 0{iJ?), ■ ■ ■ . Prom 0(/i°) the dispersion relation follows: 

uj'^ = ghk"^, or k^ (10.7) 
^/gh 
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Thus the local wave number and the local depth are related to frequency according 
to the well known dispersion relation for constant depth. As the depth decreases, the 
wavenumber increases. Hence the local phase velocity 



also decreases. 

Prom 0{ijl) we get, 

or, after multiplication by Aq, 

which implies 

hence 



c="j: = ^i (10.8) 



ax ax 



{kh\Ao\^) = (10.9) 



dx 



or. 



A;/i|Ao|^ = = constant 



^/gh\Ao\'^ — constant — ^/gh^\A^\'^ (10.10) 



Since in shallow water the group velocity equals the phase velocity, the above result 
means that the rate of energy flux is the same for all x and is consistent with the 
original assumption of unidirectional propagation. Furthermore, the local amplitude 
increases with depth as 



Aqc \ h 

This result is called Green's law. 

In summary, the leading order solution is 

1/4 /> \ 1/4 



(10.11) 



C = A^ (^^^ e^'-'-' = A^ (^^^ exp (^t k{x')dx' - lu!^ (10.12) 

Homework: Study the refraction of a train of obhquely incident waves towards a 
sloping beach. The sea depth varies slowly only in x, /i = h{x). 
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A Solution for the Green function 



Let us divide the domain into two parts, corresponding to two sides of the concentrated 
forcing. On the right side x > xq, we define G{x,xo) = G+{x,xq) which must be 
governed by : 



d'G 
dx^ 



+ k^Gj^ = 0, x > xq 



(A.l) 



To satisfy the radiation condition we take the sol,ution to beTo satisfy the radiation 
condition we take 

G+ = Ce''^"^ (A.2) 
On the left side we define Xq) = G_{x, Xq), which must be governed by 

d^G. 



dx^ 



+ k^G- = 0, X < Xq 



(A.3) 



To satisfy the radiation condition we take 
Ki X — Xq we require continuity: 

G+{xq,Xq) = G-{X),X)) 

Hence 



(A.4) 



Ce 



ikoxo 



De 



-ikoxo 



(A.5) 



(A.6) 



In addition we get another matching condition ai x — xq by integrating eq(8.8) across 
the concentrated forcing 

rxo+S f (Jp'Q ~\ rxo+S 

/ < — — ^ + k^G > dx = / 5{x — xo)dx — 1. 

Jxo-S [.dx J Jxo-S 



The left-hand side can be integrated to give the second matching codition, 

= 1 



'dG+' 




'dGJ 


dx 


xo 


dx 



xo 



implying that 



ikoCe'^^'"' + ikoDe-'^"'"' = 1 



(A.7) 



(A.8) 
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Thus 



The final solution is 



G+ = —- , x-xo>0; G_ = — — , x-xo<0. 

or more compactly, 

„iko\x—xo\ 

G{x,xo)^—— (A.IO) 
More on one-dimensional waves will be discussed later. 



B On the radiation condition 
B.l Forms of the radiation condition 

In the theory of diffraction and/or radiation of time-harmonic waves, an important 
boundary condition is the so-called radiation condition due to A. Sommerefeld : 
The disturbance due to either diffraction by a finite scatterer, or radiation from a finite 
region, can only behave as outgoing waves at infinity. 
Let us first write the unknown as 

$(x,t) = 0(x)e-^'"* (B.l) 

In a scattering (diffraction) problem we let 

= 0/ + 0' (B.2) 

where 0/ represents the incident wave (e.g., 0/ = Agt^^^), and 0' denotes the distur- 
bance (scattered wave). In a radiation problem, 0' represents the radiated wave due to 
motion in the finite region. The radiation condition can be stated mathematically in 
two equivalent forms, as follows. 

1. One-dimensional propagation: 

(P\x) Aoi3±e±*'^^ as kx iboo; (B.3) 

or, 

lim ( ^ ) ^ (B.4) 

fca;^±oo \ OX I 
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2. Two-dimensional propagation: 

(t)'{r,e,z) ^ AoB(e,z)^, as kr ^ oo; (B.5) 

y/kr 

or, 

lim Vk^ (^- ik(l)'] (B.6) 

fcr^oo \Or J 

3. Three-dimensional propagation: 

<\^{r,Q,i,)^ A^\{B,i))-—, as A;r ^ ±oo; (B.7) 

kr 

or, 



lim /cr - ^ (B.8) 

Note that the radiation condition is meant only for waves that is strictly sinusoidal in 
time from ancient past. For transient disturbances, the condition should be changed to 

$'(x,t)^0, as |x| ^ oo, t<oo. (B.9) 

Physically, if the source of the wave disturbance is in a finite region, it is not possible 
for wave energy to flow inward from infinity. Mathematically, a sinuoidal wave is the 
limit of a transient wave where the local and sinusouidal forcing begins at some time 
in the past, say t = 0. For all finite the boundary condition (B.9) holds. In order 

to derive the behavior at very large time (the quasi steady state), one can solve an 
initial-bounday value problem under the initial conditions, 

0',^ = O, t = 0, Vx (B.IO) 

and the additional boundary condition in the source region, e.g., 

^ = n • Ve-^'^*, t > 0. (B.ll) 
on 

and then take the limit of large t. The explicit solution for the IBVP is not trivial. We 
shall demonstrate for a one-dimensional example where the local energy supply grows 
slowly in comparison with the wave period. 
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B.2 Transient motion of an elastically supported string 

Consider a semi-infinte string, 

pVu-TV,:c + KV^O,, x>0, t>0 (B.12) 
subject to the boundary condition, 

V{Q,t) = Ao{et)e-''^\ t > (B.13) 

where e <^ 1 and Ao(ei) increases slowly from Aii{—oo) = to Aq{oo) — A ^constant 
with the time range ujt\ — Thus at the boundeary there is sinusoidal forcing with 

a slowly increasing amplitude. At infinity the boundary condition is 

l/(oo,t)^0, t<oo. (B.14) 

Let us use the perturbation method of mulitple scales. In vieew of (B.13), we in- 
troduce two time coordinates t and ti = et. Suggested by the elementary example of 

wave groups, wc expect that the envelope to be slowly varying in space as well. Let 
us therefore also introduce two space coordinates x,Xi = exi. We now pretend that 
V = V{x,Xi;t,ti) so that the original time derivatives become, according to the chain 
rule, 

d dV dV , 
O^V ^^+,^U5: + ,|}^U^ + 2e|^ + e^^. (B.16) 



df^ \dt dtij\dt dtj df^ dtdh dq 
Likewise the space derivatives become 

d ^dV dV 

dx dx dxi ' 

9V fd d\{dV dV\ d^V d^V ,d^V , 



dx'^ \dx dxij \dx dxij dx"^ dtdxi dxl 
We next substitute these into (B.12) 

V = Vo{x,xr,t,ti) + €Vi{x,Xi;t,ti) + €'^V{x,Xi;t,ti) + ■■■ (B.19) 

and separate different orders. At order 0(e°) we have the following perturbation equa- 
tion : 
0(£»): 
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At order 0(e), we get: 

Let us take the leading-order solution to be a sinusoidal progressive wave 

Vo^ A{xi,ti)e'^^''-'^'^ (B.22) 

where 

* = ±V^-|. (B.23) 

Let the square root be real so that there is propagation. However, we do not know the 
sign of k, hence the direction of propagation. 
At the order 0(e) Eqn. (B.21) can be written 



dt^ dx'^ \ dxi dxi 

The forcing term on the right is a homogeneous solution of the left, hence unbounded 
resonance of Vi must occur unless 

dA dA 

or 

dti pk dxi ^ ^ 
This is the condition of secularity (solvability) and is a first order hyperbolic equation. 

The general solution is 

Now we apply the boundary condition (B.13) at X\ — 0, 

A{h) = Ao{h) (B.28) 

Hence the envelope is 

For this transient problem we expect that A(oo,ii) approaches zero when x — > oo 
for finite ti, according to (B.14). Since for large xi, A{xi,ti) Ao{—xi/Tk/pu;), we 
must take the positive sign for k in view of the given behaviour of Aq. Thus 

k>0 (B.30) 
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implying that the solution is a rightward propagating wave. Thus the radiation condition 
is derived. 

Note further that 

Tk 

— = Cg (B.31) 
puj 

is the group velocity in the string, thus 

= -xi/cg) (B.32) 
or, A is a slowly varying function of x — Cgt. Thus result is true for all dispersive waves. 



